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Abstract 



A nonuniform system is considered consisting of two phases with different 
densities of particles. At each given time the distribution of the phases in 
space is chaotic: each phase filhng a set of regions with random shapes and 
locations. A chaotic diffusion process intermixes these regions by varying their 
shapes and locations in a random way. To investigate the statistical properties 
of such a system, it is exemplified by a lattice-gas model. Conditions are 
analysed when this chaotic lattice-gas model can become thermodynamically 
more stable than the usual model describing a pure one-phase system. 
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1 Introduction 



According to the space distribution of particles, there are three types of stable 
matter: uniform, periodic and quasiperiodic. The uniform matter is presented by 
gases and fluids. This can be also called plasma if it consists of charged particles. 
Periodic systems are crystals and quasiperiodic ones are incommensurate crystals 
and quasicrystals [1]. 

Amorphous sohds, or glasses, are not stable, but they are metastable. Metastable 
systems, contrary to stable ones, have finite lifetimes, although the latter can be 
quite long, as it is, e.g., for the window glass. A more detailed description of 
metastable systems can be found in reviews [2,3]. 

The dependence of system characteristics on space variables can be analysed from 
the point of view of dynamical theory by analogy with the dependence of trajectories 
on time. Basing on this analogy, amorphous solids could be called turbulent crystals 
[4]. The breaking of a regular periodicity in crystals by defects can be named the 
defect turbulence [5]. And an irregular distribution of particles, which is neither 
uniform nor periodic or quasiperiodic, is called chaotic [6]. 

As far as amorphous solids, having chaotic distribution of particles in space, 
are metastable, their consideration in time is limited by their lifetime. The state, 
existing during a finite time tj that is much longer than an observation time Tobs 
is called frozen ( Tobs t/ ) . A frozen chaotic system, after its lifetime, transforms, 
as time goes to infinity, into either uniform or periodic (quasiperiodic) system. This 
can be equivalently said that a frozen chaotic system is thermodynamically unstable. 

The question arises whether only uniform, periodic and quasiperiodic space struc- 
tures can be stable, or a spatially chaotic system also can exist an infinite time. An 
analogous question has been put forward by Ruelle [4] who has asked whether there 
can exist a frozen turbulent structure with infinite lifetime. What is known now 
about the real as well as model systems offers the negative answer to the Ruelle 
question, although this does not prove, of course, that such a situation is princip- 
ially impossible. The question that we pose here is more general: we ask whether 
there can exist a thermodynamically stable chaotic system (not necessarily frozen). 

Let us concretize this question. Suppose p{r ,t) is the distribution of particle 
density as a function of the real - space vector, r , and time, t . Assume that the 
limit 

lim p( r , t) = p(r) 

t— >oo 

exists. If the resulting state of the system with the density p(r) is stable against 
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small perturbations and has infinite lifetime, it is called thermodynamically stable. 
One can define a thermodynamic potential enjoying the property of convexity and 
being minimal for a thermodynamically stable system. The latter usually has well 
defined symmetry properties governing the behavior of the density pir) as a 
function of space coordinates [1,3,7]. If the system is uniform, then the density is 
translationally invariant with respect to an arbitrary space vector r e , that is 

p{r) — p(0). (uniformity) 

A periodic structure, or a crystal, is described by a basis {5^ | a = 1,2,3} of 
generating lattice vectors ba , with respect to which the density is periodic, 

p(r) = p(r + ba), (periodicity) 

implying that p(r) is not constant. 

A quasiperiodic structure is characterized by a basis {6a|Q; = l, 2,. 

p > 3} of integrally independent vectors b a such that the density can be expressed 
as a p - fold Fourier series 

p(7) = p„i...n^e'("i^i+'*2^2+-+"''^'')" , (quasipenodtcity) 

n\...np 

in which Ua are integers and the set {9a} is the reciprocal basis defined by the 
equations 

~gaba='2n (a = l,2,...,p> 3). 
The integral independence of ba means that the linear combination 

ni bi +n2 62 + • • • + np bp^ 

is nonzero for any nontrivial p - tuple of integers (ni, n2, . . . , Up) , i.e. when at least 
one of Ha is nonzero. Note that the structure is quasiperiodic if and only if the rank 
p of the basis {6 a} is greater than the space dimensionality: p > 3 . Although, 
probably, the majority of quasiperiodic structures are metastable, they also can 
indeed form stable equilibrium states [8] , with free energies lower than those of the 
conventional crystalline structures. Recall that if p(r) is chaotic in real space, as 
it is for different glasses, e.g. metallic glasses, vitreous silica or amorphous polymers 
[9,10], then such a system is always thermodynamically metastable. 
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Suppose now that the particle density does not tend to a stationary hmit as time 
tends to infinity, 

p(r,t)^Pc(r,t); t ^ oo, 

but displays spatiotemporal chaos [5]. The average density of the whole system is 
assumed to be fixed, 



^ J p{r,t)dr^ ^ J p^{r,t)dr^ p; 



V 

here V is the system volume. A function p{r ,t) for a given r , exhibits temporal 
chaos if its time motion is locally unstable. We may say that the motion is locally 
unstable, or locally chaotic, if a local Lyapunov exponent is positive: 



A(r ) = lim hm - In 

*^°°5p(7,to)^0 t 



6p(r,to + t) 



Sp{r,to) 



> 0. 



The case when p{r ,t) , under fixed r , tends to a time - independent limit 
p(r) , as i — > oo , describes the passage of a system to an equilibrium state. In 
the language of dynamical theory such a state corresponds to an attractor of fixed 
- point type [3]. Thus, the equilibrium uniform, spatially periodic and spatially 
quasiperiodic states are fixpoint attractors. 

The situation when a system tends, as time increases, to a chaotic state, may 
be associated with the existence of a chaotic attractor. In a system with dissipation 
the chaotic attractor is fractal, or strange [11]. As far as the equations of motion for 
p(r,t) are partial differential equations involving derivatives with respect to space 
variables as well as with respect to time, wc have to speak about spatiotemporal 
chaos. Such evolution equations correspond to what is called infinite - dimensional 
dynamical systems [12-14]. 

A chaotic state can be averaged with respect to time. For example, for the 
particle density we may define 



- ^ 1 n ^ 

p (r) = lim - / p(r,t)dt. 

r^oo r Jo 



The average state, independent of time, can be put into correspondence with an 
effective equilibrium state, or the state equilibrium on the average [3], for which 
an effective free energy can also be defined. When the latter is lower than the free 
energy of any equilibrium state, whether uniform, spatially periodic or spatially 
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quasiperiodic, then such a state is thermodynamically, or globally, stable, although 
dynamically, or locally, it is unstable. 

The question we suggest in this paper is as follows: Can a statistical system with 
spatiotemporal chaos be globally stable! We construct and analyse a model giving 
the positive answer to this question. 

2 Two - density system 

Before considering a particular realization of a chaotic statistical system it is worth 
paying some attention to terminology. Since we are dealing not with simple dy- 
namical systems but with statistical systems containing infinitely many degrees of 
freedom, we can call the chaos appearing in the latter the statistical chaos. This 
chaotic state evolves in a system automatically, without any special fine - tuning 
of external fields; hence this can be named the self - organized chaos. The time 
duration of chaotic fluctuations and their spatial sizes are distributed in a quite 
large scale, somewhat similar to microscopic critical fluctuations of equilibrium sys- 
tems, thus, one could use the term the critical chaos. This should be distinguished 
from the self - organized criticality [15] which is a generalization of the equilibrium 
critical state to nonequilibrium systems. The main feature distinguishing the chaos 
we consider is that, although each chaotic fluctuation can have any duration, size, 
and shape, but there can be deflned, nevertheless, their characteristic, or average, 
lifetime, t/ , and linear size. If , satisfying the inequalities 

Tloc <^Tf <^ Tobs, Oo < // < L, 

in which Tioc is called the local equilibrium time; Tobs , the observation time, is 
the average distance between the objects composing the system, and L is a linear 
size of the latter [3]. Because of this intermediate, or mesoscopic, character, it is 
logical to title the considered phenomenon the mesoscopic chaos. One of the most 
important observable quantities characterizing the properties of s system is the order 
parameter, generally, being a function of space and time variables. The examples 
are: the particle density p(r,t) discussed above, the local magnetization, polar- 
ization or superconducting gap. In chaotic systems, the mesoscopic fluctuations of 
order parameters are clearly interpreted as fluctuating nuclei, or gems, of competing 
phases. Therefore, such a mesoscopic chaos can be called as well the heterophase 
chaos. 
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Now, let us return to the concrete case when the chaos in a system is related to 
mesoscopic fluctuations of the particle density p{r ,t) . Suppose that at each fixed 
time the function p{r ,t) describes a random distribution in space of two regions 
with different densities, pi and p2 , such that 



where N^, and Vi, are the total number of particles and the volume, respectively, 
corresponding to a u - phase { u — 1,2 ). One can distinguish the regions occu- 
pied by different phases owing to the mesoscopic nature of the density fluctuations 
considered. The latter are accompanied, of course, by the fluctuations of other phys- 
ical quantities, for example, by rare large energy fluctuations [16,17]. However, we 
shall speak in what follows mainly about the mesoscopic density fluctuations since 
it is just density that plays the role of a natural parameter distinguishing different 
phases. 

A physical meaning which we could ascribe to these phases with two different 
densities is by calling the denser phase liquid and the less dense one gas. Then, we 
could imagine that our system is a kind of fog consisting of liquid droplets inside gas. 
Imagining this, we should not identify our chaotic fog with the real fog corresponding 
rather to a frozen droplet configuration, as the average lifetime of real droplets is of 
the order or longer than the observation time. As is discussed in the Introduction, 
such a system with a frozen space disorder is metastable, as the real fog is. However, 
some thermodynamic peculiarities of our chaotic fog can be related to the so - called 
droplet singularities [18,19] of liquid - gas systems. 

When a persistent chaotic state develops in a system, this means that such a 
state is a chaotic attractor. In dissipativc systems, chaotic attractors are fractal, 
or strange [11,12,20]. A chaotic attractor, being metrically indecomposable (met- 
rically transitive) allows to define an invariant probability measure. The property 
of metric transitivity for the considered case of a system with heterophase chaos is 
what has been called [3] the heterophase quasiergodicity. For a metrically transitive 
system, time averages can be replaced by ensemble averages. To this end, we have 
to introduce a probability measure describing the distribution of phases in space. 
The differential probability measure can be written as 



where piodO ^ locally - equilibrium density matrix, and is a functional 



Pi > P2 




V 



(1) 



^p{i) = Pioc{i)Di. 
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differential over a set 



e = {a^)k = l,2,...;7e V} 
of the manifold indicator functions 



1, rev, 

0, r^Y,; 



V and Vi, meaning the whole volume occupied by the system, and a region of 
the volume filled by a z/ - phase, respectively [3]. The averaging over the ensemble 
of phase configurations makes it possible to define an effective phase - replica sys- 
tem that is equilibrium on the average [3]. The procedure of overaging over phase 
configurations with all foundations and details has been expounded in a series of 
papers [21-24] and reviewed in [3]. As a result, we can define an effective phase - 
replica Hamiltonian which for the two - phase case takes the form 

Hi®H2, (2) 
H^^w^j (?) K{r) - ji tpt,{7)d r + 

-w^J (r) (r )$(r - r )'0i/(r )ipu{r)d r d r , 

where K{r) is the kinetic - energy operator including external fields, if any; /i , 

chemical potential; $(r — r ) , interaction potential; and '^v{f) is a field - 
operator representation for the v - phase. The renormalizing factors Wy = V^/V 
are the geometric probabilities of the corresponding phases, satisfying the condition 

1^1 W2 = 1, < < 1, (3) 

and given by the minimization of the thermodynamic potential 

y^-^XuTve-P^ (/3e = l), (4) 

in which N is the total number of particles in the system, and G is temperature. 
The total number of particles is 

f)T-f 

N^N^ + N^; 7V, = -<_^>, (5) 
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where < . . . > implies the statistical average with the phase - replica Hamiltonian 
(2). The number of particles in each of the phases is 

N^^w^J <V^^ {r)ilj^{7) >d7 . (6) 
The mean density of particles in the system is a linear combination 

p = — = wipi + W2P2 (7) 

of densities 

1 /■ t ^ ^ ^ 

Pu= y J <iljAr)Mr) > d r, (8) 

which, by definition (1), are different in the two phases. 

For a system equilibrium on average we can define all observable quantities as 
standard statistical averages with the Hamiltonian (2). Thus, for the internal energy 
per particle we have 

E=j^ (<^> +pN^ ^j^{<Hi> + <H2>)+p, 

< >= (^K, - N + wlB.N, (9) 
where the notation is used for the mean kinetic energy 

Ku = ]^j <k (r ) >d7, (10) 

and mean potential energy 

1 /■ t ^ t ^/ ^ ^/ ^ ^/ ^ 
" 27V y ^^'^ ^'^ ~ ^ )Mr)M'r )>dr (11) 

per particle in each of the phases. 

Recall that speaking here about the system equilibrium on average we imply an 
effective system obtained after averaging over phase configurations which is equiv- 
alent to the averaging with respect to time. The actual system describing chaotic 
motion of two phases is, of course, nonequilibrium [3] . 
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To define the phase probabihties, , providing the minimum for the thermo- 
dynamic potential (4), under the normahzation condition (3), it is convenient to use 
the notation 

w = Wi, W2 = 1 — w. (12) 
Then, the extremum of (4) with respect to w imphes 



^ = 1<M>=0. 

dw N dw 

This, with the Hamiltonian (2), yields 

2B2 + K2-Ki-{ii/p){p2-pi) 



w 



(13) 



(14) 



Equation (14), in compliance with notation (12) and assumption (1), defines the 
probability of the dense phase. The probability of the diluted phase is W2 = 1 — w . 
The phase probabilities given by (12) and (14) provide the minimum of (4) if 



> 0. 



(15) 



The latter inequality is equivalent to 



<^>_1< >^>o 



(16) 



Conditions (15), or (16), are conditions of heterophase stability. The second term in 
(16) is always non - negative, thence we can write a simplified version of (16), 



> > 0, 



(17) 



which is a necessary conditions of heterophase stability. The latter, taking account 
of (2), yields 

Bi + B2> 0. (18) 
In addition, the second property in (3) gives 



^ (Pl - P2) - 2Si < - ^2 < + ^ (pi - P2) 

P P 



(19) 
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Besides, it is worth checking that the system satisfies the usual conditions of 
thermodynamic stabihty, according to which the specific heat and isothermal com- 
pressibily are to be positive: 

-1 



= [p^ j > 0, (20) 

here the free energy per particle, / , and the pressure, P , are related to the 
thermodynamical potential (4) by the equations 

f = ye + P, 

dp 

The system is called absolutely stable in the thermodynamic sense if its thermody- 
namic potential is the lowest possible and all stability conditions described above 
are valid. For brevity, we may call this type of absolute thermodynamic stability the 
global stability, as compared with the local dynamic instability of chaotic motion. 
The heterophase chaotic system is dynamically (locally) unstable, while it can be 
thermodynamically (globally) stable. 

3 Lattice - gas model 

To get more details of statistical properties for a system with heterophase chaos, 
let us pass to an approximation leading to what is termed lattice - gas model. To 
this end, we divide the whole volume of the system, V , into Nl lattice cells, the 
density of lattice sites being 

Po ^ (21) 
This quantity can be used for defining the dimensionless densities 

_ Pu _ P N 

nu = —, n= — = —. 22 
Po Po Nl 

With notation (22), we can recast (7) to the equation 

n = Wifii + W2n2, (23) 
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which may be used for defining the chemical potential as a function of density 
and temperature, /i(n, 0) . 

Introduce a set {a, \ i = 1,2, . . . , N^} of vectors corresponding to lattice sites 
and define the Wannier functions (pnif — ^i) with n indexing quantum states. 
The field operator ipu{f) can be written as an expansion 

'^A^) = Cinueiu^n(r - ai), (24) 

i=l n 

in which Cji, = 1 , is a variable telling whether the i - lattice cell inside the - 
phase is occupied by a particle or empty, respectively. Assume that each cell can 
accept not more than one particle, which results in the homeopolarity condition 

n 

A standard supposition for a lattice - gas approximation [25] is that intercell and 
interlevel transitions are prohibited. This can be expressed in two ways. One pos- 
sibility [26] is to require the diagonality of the matrix elements 

<mi\K\nj >= 6„inSijKi, 

< mi,m'i'\^n'j' ,nj >= SmnSijSm'n'Si'j'^ij. 

Equivalently, we can define the operators Cinu on a restricted space of states with 
the restriction 

Then, (24), (25) and (26) permit to cast each term of (2) into the form 

H„ = - i^Yiy + ^ij^iuejy. (27) 

Thus, we obtain the effective phase - replica Hamiltonian H= ifi © , with terms 
(27), describing the statistical properties of a chaotic lattice - gas model. 

In the considered case of heterophase chaos, the system consists of two randomly 
distributed phases with different densities. The latter are given by (8) which for the 
lattice gas become 

I Nl 

= T7 E < eii. > ■ (28) 

^ i=l 
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As the dimensionless densities (22) we get 



I Nl 



1=1 



Nl 



E < > • (29) 



From the occupation - of - cell variable Cj,^ = 1, one can pass to the spin variable 
— ^1 by using the transformation 

eiu ^ ^ (1 + ^i'^) ' ^i'^ ^ 2ei^ - 1. (30) 



Then, for the density (29) we have 



n^ = ^{l + mi,), (31) 

where 

I Nl 

^ < (Jiu > ■ (32) 

According to (1), the first phase is dense and the second one is diluted, which for 
(31) and (32) means that 

rii > 712-1 TTli > 777.2. (33) 

In the quasispin representation the Hamiltonian (27) becomes 

I Nl I Nl 

with the notation 

I Nl I Nl 

Nr ^ " Nr ^ 

hiy = n - Ki - 

The operator from of (34) corresponds to the Ising model with an external field. 
Such a model, as is known, is difficult to treat. For low dimensionality, there exists 
a numerical transfer - matrix method based on a real - space Trotter decomposition 
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of the statistical operator [27]. A numerical transfer - matrix study based on finite 
- size scaling [28] is applicable to the Ising model in two and three dimensions. 
However, our aim here is not to present some extra accurate calculations but to 
clarify the main physical features of a new model. Therefore, at this stage it would 
be unreasonable to sink into clumsy technical details. Instead, we can think that 
$jj is a long - range interaction for which the mean - field approximation works 
well. The latter gives simple expressions for correlation functions, such as 

<ai„aj^>^ml (i^j)- 

For arbitrary lattice sites we can write 

< (Ti^cTj^ >= Sij + (1 - 5ij)ml. (35) 

Using (35), we are able to continue the explicit analytical investigation of our model. 

First of all, let us check the conditions of heterophase stability. For quantities 
(10) and (11) we now have 

K , , n,, 

Jf„ = _(l + ,„„) = A-^. 

B,= *(l + ,„„f = *|. (36) 
As far as in the mean - field approximation 

2 

> = < 0, 

aw 

the conditions of stability (16) and (17) coincide with each other and with (18). The 
latter, taking account of (36), leads to the condition of heterophase stabihty 

$ > 0. (37) 

This shows that to make the heterophase system stable, the particle interactions are 
to be repulsive on average. 

For the free energy we obtain 
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n V 



2 cosh 



$(1 + m^)wl + 2{K - ii)wy 



40 

The mean quasispin (32) becomes 



(38) 



, ^ ^wl{l + m,) + 2w,{K-ix) 

mi, = — tanh -— . 39 

4© 

The equation for the phase probabihties can be obtained either from(13), that is 
from the minimization of (38) with respect to -u; , or by substituting (36) into (14). 
In the both cases we get, of course, the same equation 

wr^{l + mif + 2mi{K - //) = W2^{1 + ms)^ + 2m2{K - /i), (40) 

to which we have to add (23) to define the chemical potentiaL 

4 Thermodynamic properties 

Now we have in hands all information needed for investigating thermodynamic prop- 
erties of the chaotic lattice - gas model, whose phase probabilities and chemical po- 
tential are given by equations (40) and (23). For convenience, we write down these 
equation s in the form 

_ ^nl + (K-ii)(n2-ni) 

n = wni + (1 — 'w)n2- (42) 

First, let us show that the considered chaotic state is thermodynamically more 
profitable than a homophase state. For zero temperature, this can be readily shown 
by comparing the internal energies of the chaotic two - phase system and of a 
homogeneous one - phase system. Thus, substituting (36) into (9), we find 

K $/22 2 2\ 

£/ = — (wirii -\- W2n2) + — [WiTii + W2n2j ■ 
Invoking (42), we obtain 

n $ 

E ^ K + $ wiW2nin2. (43) 

2 n 
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In the case of a homophase system, when one of the phase probabihties is zero, we 
have 

E,,,^^K + ^-. (44) 
The comparison of (43) and (44) proves that 

E < Enom (^ > 0). (45) 
At finite temperatures, we have to compare (38) with 

= + + In (2cosh ^) , 

where 

K + (^n + 2Qarctanh{2n - 1). 
Numerical analysis shows that there always exists a range of temperatures when 

f fhom ■ 

To investigate in more details the temperature behavior of thermodynamic char- 
acteristics of the chaotic lattice - gas model, and to check its stability conditions, 
we proceed in what follows to a numerical investigation. For this purpose it is 
convenient introduce a dimensionless chemical potential and temperature, 

T.f. (46) 

Let us also use the notation 

a = rii, b = n2 {a > b). (47) 
Eq.(38), with the help of (46) and (47), can be rewritten as 

= ^ (1 — a) -u;^ + ^ (1 — 6) (1 — w)^ — -fji*— 



$ 2 ' ' 2 ' " '2 

Tin 



, aw"^ — a^w , bil — wY — LiAl — w) 
4 cosh — f- cosh ■ / V 



2T 2T 



(48) 



The quantities a and b play the role of order parameters distinguishing phases 
with different densities. The equations for these order parameters follow from (39) 
giving 

2a-l^-tanh"^'-^-^ 
2T 
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26 - 1 = - tanh 



6(1 - w)^ - - w) 
2T ' 



(49) 



And the equations for w and are (41) and (42). From (41) we can express the 
chemical potential 

a^w — 6^(1 — w) 
A** = - 



a — b 



and from (42), the phase probability 

n — b 
a — b 



w 



(a > b). 



(50) 



(51) 



Substituting (51) into (50), we get 

a^(n — 6) + 6^(n — a) 
(^^^6p ■ 

Equations (49), with the use of (52), become 

2ab — n{a + b) 



(52) 



2a = 1 — tanh 



26 = 1 + tanh 



b{n - b) 
a{n — a) 



2T{a - 6)3 
2a6 — n(a + 6) 



(53) 



2T(a - 6)3 

In this way, if we solve (53) with respect to a > b , then we can calculate (51) and 
(52). 

Note that (51), together with the property of the phase probability, < ty < 1 , 
impose the following condition 

b < n< a (54) 

on the solution of (53). Numerical investigation of (53) shows that the solutions 
a > b exist provided that 

< n< 0.323. (55) 

These solutions are plotted in fig.l. 

Figs. 2 and 3 display the temperature behaviour of the dense - phase probability 
(51) and of the dimensionless chemical potential (52), respectively. As we see, the 
probability of the dense phase diminishes with temperature, becoming zero at some 
Tc{n) . 
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The pressure of the system is drawn in fig. 4. And figs. 5 and 6 demonstrate that 
the stabihty conditions (20) for the specific heat and isothermal compressibihty are 
satisfied. 

To illustrate that the condition of heterophase stability (37) is directly related to 
the conditions of thermodynamic stability (20), we also considered the case $ < . 
The corresponding equations for the order parameters can be obtained from (53) 
after the change $ — >■ — $ . The resulting equations do possess nontrivial solutions 
depicted in figs. 7, 8 and 9. Even more, there appear three branches of these solutions. 
However, for no one of these branches the conditions of thermodynamic stability (20) 
hold. Thus, the conditions of heterophase and thermodynamic stability are either 
both valid or both broken. 



5 Discussion 

We have constructed a heterophase lattice - gas model consisting of two phases with 
different densities, one phase may be called dense; another, diluted. The dynamics 
of the phases is chaotic. However, the system, as a whole, can be globally stable. 

The density of particles in such a system, after averaging over phase configura- 
tions, or over time, has the form of a sum 

^ (^) =E^-P-(^)' (56) 

in which 

^ t ^ ^ 
p^(r) =<^^ (r)V^^(r) > . (57) 

The space averages of (56) and (57), i.e., 

lead to the mean densities (7) and (8), respectively. 

For the lattice - gas model considered, the density (56), in accordance with (57) 
and (24), is periodic. In this way, although at each phase configuration, or at any 
fixed time, there is no periodicity in the system because of random distribution of 
two phases is space, but after averaging over phase configurations, or over time, the 
averaged density of particles (56) happens to be periodic. 

A straightforward generalization of the model could be if we would ascribe to 
different phases different lattices. Putting into correspondence to each phase Nli, 
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cells with the attached lattice vectors {ai^} , we may write, instead of (24), the 
expansion 

= Y,^iriueiuVn(r - a). (59) 

i=l n 

Then, in place of (57) we would have 

^ t 
PA^) = < ^i'^ XCinu Cinu > \^n(r - a^i.)|^ (60) 

i=l n 

which is a periodic function, 

p^(r + ai^) = p^(r) 

If for a two - phase system the vectors Oij, and are integrally independent, 

as is defined in the Introduction, then the sum (56) of two periodic functions (59) 
is a quasipcriodic function. Similarly, wc could construct a quasiperiodic linear 
combination (56) composed of an arbitrary number of periodic terms. 

So, a system can seem, on average, to be quasiperiodic while actually it is chaotic. 
It is possible to speculate that some incommensurate structures could be of this 
kind. Whether an order in a system is frozen or is sustained only on average, can 
be checked by appropriate experiments. For example, the Mossbauer spectroscopy 
is very sensitive to structural features, and the Mossbauer factor displays specific 
anomalies in the presence of lattice instabilities, especially when the latter are con- 
nected with mesoscopic structural fluctuations as in the considered case [29,30]. 
Anyway, even if such systems would not exist in nature, the model investigated here 
presents an illustration of unexpected possibility when the onset of chaos, that is 
of local instability, makes the system globally more stable. However, we think that 
this possibility is not just a strangeness of an artificial model but rather a general 
phenomenon common for many physical systems, the phenomenon called [3] the 
spontaneous breaking of equilibrium. 
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Figure captions 



Fig.l. 

The dimensionless density of the dense phase (sohd curves) and the diluted phase 
(dashed curves) vs. temperature for several mean densities. The curves end at the 
points where the system looses its stability. 

Fig.2 

The probability of the dense phase as a function of temperature. 
Fig.3 

The dimensionless chemical potential vs. temperature. 
Fig.4 

The temperature dependence of pressure. 
Fig.5 

The specific heat as a function of temperature. 



Fig.6 

The isotermal compressibility vs. temperature. 
Fig.7 

The dimensionless density of the dense phase (sohd curves) and of the diluted phase 
(dashed curves) vs. temperature in the case of the unstable system with the mean 
density n — 0.20 . Numbers 1,2 and 3 enumerate the branches of solutions. 

Fig.8 
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The same as in fig. 7 but for tlie mean density n = 0.26 . 
Fig.9 

Tlie dependence of tlie densities of the dense phase (sohd curve) and of the diluted 
phase (dashed curve) on the mean density in the case of the unstable system at fixed 
temperature T = 0.03 . 
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